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ABSTRACT 


The general equations of surface ship motion are 
developed and standardized for simulation in digital 
computer. Digital simulations of the dynamics of the 
surface ship in three degrees of freedom are done 


with and without non-linear and cross-coupling terms. 
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ก โก พ ร ค ก κενο SHIP MOTION 


noa Ship ie E MEbESIx degrees of freedom. 
These degrees of freedom are generally chosen as follows: 
near απο πο along the three axes 


merough the center of gravity. 


a.l Surge -- along X axes 
a.2 Sway -- along Y axes 
a.3 Heave -- along Z axes 


b. Rotations around the three axes through the 


Ecnter of gravity. 


b.l Roll  -- around X axes 
b.2 Pitch -- around Y axes 
b.3 Yaw -- around Z axes 


Further reduction in the complex nature of the equa- 
tions can be brought about by choosing an orthogonal 
axis system parallel to the principal axes of inertia so 
as to eliminate products of inertia in the motion equa- 
tions. For practically all ocean vehicles, with extreme- 
ly few exceptions, a longitudinal axis (X axis) in the 
centerline plane, a downward (toward keel) axis (Z axis) 
perpendicular to the X axis in the centerline plane, 
and a transverse axis (Y axis) perpendicular to the 


centerline plane satisfies this requirement. 








+ HEAVE 


Figure l. Surface ship in six degrees of freedom. 


For the exceptional vehicle which has a very peculiar 
and significantly large asymmetrical mass distribution, 
it is necessary to include the products of inertia. 

КОО c E 7 axes form an orthogonal right hand system 
of axes fixed in the vehicle. The axes and the associated 
components of the pertinent physical quantities are 
defined below: 

The longitudinal X axis (in the plane of symmetry) is 
positive in the forward direction, usually parallel to 
the keel or calm water line. If the upper and lower halves 
of the body are symmetrical, then the axis is the inter- 
section of the two planes of symmetry. 

The Y axis is the transverse axis perpendicular to the 
plane of symmetry and positive to the starboard. 

The Z axis or downward axis in the plane of symmetry 
(X,2) is perpendicular to the X axis and positive downward 


towards the keel. 





Mm) nit vectors along the X,Y,2 axis respectively. 
R = tr distance Of a point from the origin O, 
and the corresponding components along the X,Y and Z axes. 


R = i Fg + k 
ου ЕС 


б u,v,w velocity or the origin O (on the body) and 


the corresponding components along the X,Y and Z axis. 


U = iu + jv + kw 
ὦ p,q,r angular velocity of the body about the origin 


and the corresponding components about the axes. 


^ ^ ^ 


б = า ว + jn + k 
The moments of inertia of the body about the X,Y and 


Z axes respectively І, Lon І,. 


F X,Y,Z, force acting on the body and the correspond- 


ing components along the axes. 


—y 


F = iX + jY + kZ 


m K,M,N Moments acting about the axes. 
> 2 A - 
m = iK + jM + KN 
Newton's law of motion for a rigid body can be written 
as two equations, one a force equation and the second a 
moment equation provided an origin is taken at the center 
of gravity and the axis system is fixed in space. The 


equations are 
> 


>m а κα > 

F = dt (momentum) T (m Ug) 

m = Sie, ular CM. = d (18) 
dt 9 G dt 








where the subscript G refers to an origin at the center 
Of gravity and m is the mass of the body. For a mass 


essentially constant in time 


> > 
F = m (UL) 


E 
ct 


For an origin not at the center of gravity of the 
body and in a system of axes fixed in and moving with 


Bre vehicle. 


where U. is the velocity of the origin in space. However, 
Since the origin is on the surface of the earth and the 


earth rotates, then 


where U is the geographical velocity of the body, ü ls 
the angular velocity of the earth, and E is the radius 
vector from earth's center to the vehicle. The force 


equation becomes: 


— d > > > 
F = m gg (U + à me Rp + Q x R.) 
> | а ,> > > > > + 
F = № qg (U + Q x Ry) +m[ e R +Q x В) 
E а ,> > > 
= 0 ag (U + ὦ x Ас) но (б ке +0 x NE Rod 


since 5 - 0 апа Ro 


ll 
Civ 
+ 

1 
ж 
0 


> > . : 3 -> 
the term ณ์ 2 ХО τηε ๑ ร ท า ๑ 115 force апа mb x ü X Ry 


is the centripetal acceleration due to rotation of the 





earth. These two terms are negligibly small when com- 


pared with the other forces, then 


О, 


> > 
Е = тю ар (9 κ Ὁ κ Re) 


Finding the derivatives of unit vectors (change in 


direction) 
where di = -kd9 di = jal ді з 0 
dj = 0 dj = -id dj = каф 
dk - id dk = 0 dk = -јаф 
Adding the contributions 
πα. ^ар 246 
ac ^ 17390 9 Зак 7 Кае 
Cu ου... , d$ 
dt ^ iac * 230 + Kae 
SINE do NM M ^ 
dt ^ Ll dc ^ Jac * F^? 
Q = ip + jq + kr 
_ аф GO _ dv 
Pr I!" 
di =S 3 ет τ 2 k 
HIS X 1 = р q ү 
1 0 0 
= ο `. 5 
ш = E ку 
1 EL 
ο οσα να с 
au А ^d: А 42 A dk 
к = ji = y a х = 
dt πω. “Ост ΠΤ. “ас 
_ te ου зае “ау , “аф 
1 น + «бак kay) НЕГІ! It και) 
as ^dO “аф 
+ Кы + м (і Jae 


= iu+ u(jr- kq)+ jv+ v(-ir+ kp)+kw + w(ig - jp) 
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i j K 


> а a τ 
Q x R. = | p q r - 7 ร 7 ก YX IF (PY ахо) 
^e Yema 
а = ENS s = E 2 Ф E 2 
Te (ÜxR = і (24 Yar) + J (2 qr Y ¿Y ) + k (Y ar 2-4 ) 


ñ . . τ 2 Š 
r j (XoY = 26р) + к (Хогр 26р ) + і (2 Tp Mp) 


ο ‘ : À 2 Е | 
I M ไธ โร”+”"”"+ - хер 
from F = ix + JY + kZ and É = WES * ac x R 
Rewriting and grouping all i terms equal to X, all j 


terms equal to Y and all k terms equal to Z yield 


Е - : D e 
X = m(u vr+wgt2q + rpZc r Xot PAY o q Xo? 
. o 
= mM (u+wg-vr=X¿ (1 за") + ς (ρᾳ-τ) + Z. (rp+q) 
° ° 2 2 ° 
Y = m (vtur-wptXe (pq+r) - Ус (Р зх") + Z5 (pa-r) (1) 
. . 9 2 2 
= m (w*vp-uq*X. (rp-q) + Y. (qr+p) - Za (P +q) 
F m, = E 1 | 
rom m. = artangular momentum) e 
on 


O, 


a+ ki, r) 
t Xe Үс Za 


- ๑ ๕ am Origin at the center of gravity 


= a 2 2 
I, = І, m Yo + 26 ) 
е 
_ 2 2 
I. = L, m (2 + Xo ) 
G 
- 2 2 
I, = I, m (Χς + Ус ) 
G 
> > > > > + > -> 
M = Me + Re x F ος Ме = M- Re ΧΕ 


ΠΠ 





After manipulating and using the results for the 
derivatives of unit vectors (same as above) expressions 


for K, H and N are obtained. 


K = І Р +(I,- I, )qr+m[Y. (w+ pv - qu)- Zg(v*ru-pw)] 
M = Iq t, - I,)rptml2,(u- qw - rv)- X¿(Ww+pv-qu) ] (2) 
N — I,r *(L,- I.) pq+m[X. (v+ ru - pw)- Yg (utqw-rv)] 


The terms (qw-rv), (ru-pw) and (pv-qu) are gyroscopic 
effects. 


ιν l o on bip i orrrEorces and moments can be ex- 


pressed 
Force) 
) = £(properties of body, properties of 
Moment) motion, properties of fluid) 


Бо Ка Бат сп а сїр, їп а given fluid with no 


excitation force - so 
Fo) 
: า ts of fetion) 
М ) 
- £(X,,Y,,2,,0,0 ju, V, Ww,p,Q,r,u,V,w, 


p,q,r,6,6',6") 
The Taylor series which has the following form may 
now be applied to linearize the equations about an 


operating point Хо 
อ อ อ อ ล ล ว ล ฑ ๓ )+ 
๐ о ах O O сьо з 


Арріу this бо аы 


12 





For the case of £(properties of motion) let (XX) 


= АХ, Nu ve) = AY. and (2-2) = А2 


terms second order and higher are neglected for small 
perturbations. 


From X equation the linear terms are obtained: 
- + AX. EE xU AY 25%. ze e 8 : 0 m 


The defining relations are: 


«Ое 7 ‘Sa? usu, en 
6X " 
(วั น พ = พ З 
Aw = (w-w_)=w, vu - 0 
๐ ๐ 
Au = (στα) 


The force equations then become 


X = Хо ШЫ YY Mere, ด 0 200 A 


Хо O Хо O 0 u u 
Y = YON Ша Хо + ne Y %....Ү40 + eee Yu (3) 
๐ : Ο 
227 Í Zx Xo E ty Yo +.… Z8 + e AA 


A similar derivation can be done for the K,M and N equa- 
tions. The preceding X,Y and Z equations may now be 
equated to the linearized X,Y and Z (equations (1)), e.g., 
for the Y equation without roll, pitch, and the center 
of gravity at Ха 07 Үст 0 and ZT 0 gives the linearizeq 


equation. 


Y = m(v + ur) 
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then 


Yy Í Yy τοῦ Y 
O O 


Y+ Y ut Y vt Y r+Y-u+ Yevt Yer = m(v+ru) 
u v qe u V т 


ф 


Expressions for X,Z and K,M and N can be determined in 
a Similar procedure. 

In order to obtain the equations in a non-dimensional 
form some definitions will be given, and applied to the Y 


force equations as an example of the process. 





Froude number = a 
vg! 
u Уг 
/4/3 
bu AA 
X Y, 
λα αμ ως ว้ 
042 
! ! ' = ก ก 
Ea c 1 2 2 


After replacing and adding the effect of waves, the 


Y equation becomes: 


vit r'u! = 5 0 V αρ Cy F т Cy r'+ С ถา ) 
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πι КЕ ОТЕ SIMULATION 


The six equations of motion after rearranging by plac- 
ing the second order terms to the left and the rest on the 
right side become: 


aal+baB+caC+daD+eaE+faF = -(a 1a] Ataya,Atbja]Btb, a 28 


+ срауС*сзазс+ајаур+4 азр 


eja ,Ete,a, E+f a F+f a F) 


+ โท 1 


十 


abA+bbB+cbC+dbD+ebE+fbF = -(аувулна, b ΩΣ, 


+ yb, Cte yb C Hd b, Dé bn 


+ e,byEte,b_E+£,b,P+£,b, F) 


TPIP 2 


acA+bcBt+ccC+dcDtecEt+fcF = -(аусулча, c Mb, e Bib, 


+ 059405656, e P«d,c,D 


ісі Еғе, e jE*f c. F+£ >C 


+ IF3 


+ e F) 


2 


ума, а jb, di Bib йв 


+ орбустоздусна дурна, ар 


+ e,d,Ete,d,E+f, d, F+£,d,F) 


+ IF4 


adA+bdB+cdC+d@D+edE+faF = - (a, d 


аеА+реВ+сеС+аер+ееЕ+ҒеЕ = Ata, e A+b.e B+b e B 


- (ae, ภพ 


+ c,e,C+c,e,C+d, e,D+d,e,D 


12, E+e,0,E+f,e,F+f,e PF Р) 


EIFS 


+ е 
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afA+bfB+cfC+dfD+efE+ffF = -(a,f,Ata,f,Atb, f)B+rb,f,B 


мога c,f£,C+c,f,ctd, f,D+d,f,D 


+ e,f Ete f,Etf,f F+f,f,F) 


355 10046 
where A=u, A=u, Ev В-у, Сем, Сем, бер, рер 
E=q, Е-а, F=r, F=r, terms IF include all non-linear 


terms such as wave force, wind, rudder deflection---etc. 


In the six equations, the right can be set equal to 


Τι, І., ------ Іс respectivelv, thus: 

โว = -(a,a,Ata,a,Atb,a,B+ 一 一 一 一 一 一 一 一 一 一 f,a,F) + IFl 
I, = -(a,b,Atajb,Atb,b, Bt ---------- £,b,F) + ТЕ? 
โว = - (a, c,Ata,c,Atb, C, B+ 一 一 一 一 一 一 一 一 一 一 £,c,F) + IF3 
I, = - (893 A+a d Atb} d B+ ---------- f d F) + IF4 
I. = -(a,e,Ata,e,Atb,e,B+ 一 一 一 一 一 一 一 一 一 一 f,e,F) + IF5 
I; = - (a f Ata f Atb) f} B+ ---------- f,f,F) + IF6 


the equations then have the form that follows, 


aah + baB + саС + дар + eaE + faF = I 


1 
abA + bbB + cbC + dbD + ebE + fbF = Т 
асА + bcB + cec + аср + ecE + fcF = โว 
adA + Ьав + сас + аар + edE + fdF = un 
aeA + beB + ceC + deD + eeE + feF = I. 
afA ^ bfB + cfC + dfD + efE + ffF = τ, 


16 





expressing in matrix form: 


Fh 
Ф 
>: 
H 


aa ba ca da ea 


її 
ab ЬЬ cb аЬ eb £b В ix 
ac BE ee ас ес re Ç โว 
ad ва са dd ed ға ก т, 
ae be ce de ee fe Е โซ 
аї bf cf df ef fa F ы 
Apply Cramer's rule to solve for A, 1  . F in 
terms of Ii == Ic 
Ii ba ca da ea fa 
I, bb cb ар ер fb 
โว Бе C C ас ес £ C 
І, bd са ад ед ға 
โธ ре се ае ЕЕ Ғе 
Й Ic рі er dr ef ΓΕ 
аа ра са аа еа fa 


ab bb cb ар ер fb 

ac bc Ge ас ес fe 

ad bd cd dd ed fd 

ae be ce de ee fe 

аї bf GE ағ (208 ЕЕ 
define the denominator determinant 2 Жие тог the 
nominator let cofactor of Ii = col du 
Seractor of โว Γ cof. ab, --------- and coefactor of I 


A 
= cof. af equations becomes: 


6 


(cof.aal,+cof.abI,+cof.acI,tcof.adI,+cof.ael,tcof.afl,) 
A 


ju 





ШОО е сае мау solve for B, С, ---------- 


E (cof bal, СОГ ไว ไว โว cof bel, СОЕ bal, Cor bel; Cor bfl.) 

B = A 

В (сої сат, СОЁ cbI, СОЕ ссІ. СОЁ сат, cof cel. COL cfl.) 

С з น 11 — —P 

a (cof даї, COE арі, сої ๑ ๕ ๐ 1 ว СОГ дат, СОГ del. COE AfI¿) 
7 Δ 

m (СОГ eal, Cor ері. COF ес, COE еді; СОЕ eel, COF егіс) 
4 เณ ไ ไ โฆ ซ้ ซู 

„Ши к жоу л у = cof rer, eof FETS) 
q À 

meen the value of A, A, B, B --------- F, F by integration 


sien that 


A block diagram 
set of equations 15 


In the computer 


2 , 


to compute all of the variables in the 
presented in Fig. 2. 


program that is used for simulation 


all six equations for six degrees of freedom are used, but 


are interested in less than six degrees of freedom. The 


Same program can be 


used by setting the coupling terms of 


non-used equations equal to zero and one in terms of princi- 


pal diagonal, e.g. only three degrees of freedom are used 


in this study, surge, sway and yaw, then all coupling terms 


18 
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Fig. 2 


TRANSFER FUNCTION BLOOK DIAGRAM SIX DEGREES OF FREEDOM(Eqn.A) 
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Ріє. а 
TRANSEER FUNCTION BLOCK DIAGRAM SIX DEGREES OF FREEDOM(Eqn.B) 
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| 


A— 446; 
oo Ыс 
B— D.C; 





Рід. 2 
TRANSFER FUNCTION BLOCK DIAGRAM SIX DEGREES OF FREEDON(Eqn.C) 
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db 
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Fig. 2 
TRANSFER FUNCTION BLOCK DIAGRAM SIX DEGREE OF FREEDON (Eqn. D) 
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Fig. 2 
TRANSFER FUNCTION BLOCK DIAGRAN SIX DEGREES OF FREEDOM(Eqn.E) 


23 








Fie. 2 


TRANSFER FUNCTION BLOCK DIAGRAW SIX DEGREES OF FREEDOM (Eqn.F) 
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are set equal to zero and unused terms on the principal 


diagonal equal to one, for example : 


aa ba 0 0 0 fa A I] 
ab bb 0 0 0 fb B i 
0 0 1 0 0 0 С E - 
0 0 0 1 0 0 Bl Т, 
0 0 0 0 i 0 E Т. 
аї bf 0 0 0 ff F το 


and the left side of the unused equations are set equal 
Bo zero. 

With this program non-linear terms can be added such 
as, rudder deflection of waves and wind, etc., which will 
be done by adding all of these whose sum is IF, e.g. 

1F, = KA, X Dr+KA, x рѕ+КА.х Db+NA 

where Dr, Ds and Db are rudder deflection, canard de- 
MON... e. etc. NA is the sum of all non-linear 
terms that effect the surge equation (X equation). 

The program that will be used for solving these equa- 
tions is the "Continuous Systems Modeling Program" (CSMP) 


[Ref. 3) in which all constants are declared in the first 


section and then set the value of matrix for aa, ab, ac 


99» and so on (in program AAA is used for aa, AAB for 
Шен... AFF for ff). In the initial section values of 
the COFACTORS aa, ba....are determined. All of the CO- 


FACTORS and the subprogram VALUE is used to compute. This 


Subprogram finds the value of the determinant of the 


25 





matrix. For all of the COFACTOR terms the element is set 
equal to one and the rest of the elements in that row and 
column are set equal to zero. For example, to find the 
value of COF.aa the following array is obtained: 
1 0 0 0 0 0 
0 Ep co db ер” тр 
COF aa = 0 "ΕΕ ec “EC 
0 ра са dd ed fd 
0 be ce de ee fe 
0 ΕΙ dh er “тт 
(In the computer program BAA is used for 8181 САА for 
а,а,...). After the value of A and all cofactors are 
determined, the dynamic section is used to determine BAA, 
ЕСЕ ...... GAA, GAB ...... (if those terms contain variables). 
In the dynamic section all variables that are functions 
of time are determined. The defining relations of the 
variables are also included in the dynamics section, i.e. 
UDOT - ADDOT аон U=ADOT.......etc. XH,YH,ZH аге deter- 
mined and are the vector terms, Х,Ү,2 whose origin is fixed 


on the earth (relative to the earth). 
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IUE ОЕ СТР Тр PERFORMANCE 


In this section the computer will be used to solve 
the equation of motion describing ship "D". The hydro- 
dynamic coefficients and constants that were obtained 
from NSRDC (NAVAL SHIP RESEARCH AND DEVELOPMENT CENTER) 
[Ref. 6] for this study will concern only three degrees 
of freedom such that SURGE, SWAY, and YAW. 

Equations of motion of the ship (nondimensional) 


Axial Force 


m(u-vr+wg) = = Lt cg" + € + Xp YB) 
+ 5 (X. + X, vr + Хо) 
+ SA, N 2) 
ES u^ (X, Or HK, ล 8) 
+ 5 (Ауу + Ajuru + Ben 3 
Gateral Force 
m(v-ur-wp) - 5 L (Yap t Yer + ο 
+ = tp 
+ 5 (uY v 1 WV 4 ร уу) 
& By u^ (Y, Sr) 
Yawing Moment 
NI (І зі дра - 5 (Nr + Νορ + N, pa) 
+ 5; V+uN p+uN, r+N ES 
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+ Ë (uN v+N _wv+N |м |у) 
v ΨΥ 


2 


25 |м |у 


р 2 
+ —— u (N. ór) 
292 ST 


P, the density of fluid is taken as 2 and the terms 
including w, p. q (heave, roll and pitch) are set equal 
Б д2его. 

Тһе nonlinear terms such as the squared terms and 
product terms of v and r are omitted initially. This is 
in agreement with the small perturbation theory. 


After rearranging, the equations become 


: 2 u* 
Xy ล น 
2 
я 4 ° , = = u -— - uU ` 
(Y. ου Ὁ Yer τ ἴον uy r 7 Y spòr 
. 2 
М.у ο... | au _ u 
Ут ar (N.-I,) r = E — N r πο Ns Or 
x 2 2 
Set the left side of the equations equal to 1η ετττῖς, 
then the matrix equation becomes 
(X. m) 0 0. 0 0 0 А Ii 
0 (Υ -m) 0 0 0 Y, B I, 
0 0 I 0 0 0 C ІЗ 
0 0 0 1 0 0 р І, 
0 0 0 0 1 0 Е I, 
0 N./1 0 0 0 (N.-I,) F Іс 
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Also set the right side of the equations equal to I 


N 


ย ง 


> 


л 


о 


о 


о 


о 


о 


0 0 
0 0 
V 
0 0 
0 0 
0 0 
x 0 0 
Ail Zero 
BT = F, and 
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ΠΠ 
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tj 
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А. NONLINEAR TERMS NOT INCLUDED 


TABLE SI 


fe oreaynamic Coefficients and Constants of Ship "D" 
for Linear Terms (Non-dimensional) 


m = 0,0045 
Iz = 0.0003 
Nr = 0.0012 
Nr = -0.0002 
Nv = -0.0012 
Nv = -0.0001 
Xu = -0.00036 
Yv = -0.0025 
Yv = -0.0063 
Yr = 0.004 
Xsrér = —0.0011 
Үсү = 0.0019 
Nor = -0.00084 К 


6s and ób equal zero 
In the computer program, set all coefficients in 
section 1 and set AAA = Xi-m, AAB = 0 ---- AFF = Nr-Iz. 


After this, use subprogram value to find the determinant 


and coefficient of AA, BB ---- and then set BBB = T Yv , 
BFB - uYr ------- BFF = Ч Nr : 
: 2 
IFl = KAlór where KAl з -Х Gis u 
eror $, 
IF2 = KBlór where KB2 = -Үшу Шә 
IF3 = KFlör where KFl = Νες 42/42 
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Following is the program that used CSMP to determine 
the turning circles for rudder angles of SO A rad.), 
259(-0.4365 rad.) and 35?(-0.6111 rad.). Result of this 
study is presented in Fig. 3. The turning rate as a func- 
tion of time is interested and the results of this analysis 
are presented in Fig. 4. The ships turning performance 
expressed in transfer ship lengths as time is shown in 
Fig. 5 and the heading angle as a function of time is 
given in Fig. 6. Fig. 7 shows the results of the zig-zag 
maneuver, curve shown the yaw angle and rudder angle in 
degree as a function of time, for this study the same pro- 


gram was used, but set DR in dynamic section: 


DR= -0.06984 (RAMP (0.0) -RAMP (5.0) )+0.04656(RAMP(40.0)... 
-КАМР (65.0) -ВАМР (145.0) +ВАМР (170.0) +РАМР (250.0) .... 


-RAMP (265.0) -RAMP (345.0) *RAMP (425.0) * RAMP (440.0)) 


and use prepare statement prepare X, YAWD and prepare X, 


DOO (YAWD = YAW* 57.273, DOO - DR* 57.273). 
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В. INCLUDED NONLINEAR TERMS 
The computer program 2 is the same as the computer 
ΠΠ, but nonlinear terms are added by setting: 


NA = NAl + NA2 + NA3 + NA4 


where NAl = (Xx + EE + X тр) 
NA2 = -(mvr + Хх „УК + 0 + mwg) 
МАЗ = -(Х v^) /1 
NA4 = Ns + A2uu, + 0 
NB = NBl + NB2 + NB3 

where NBl = ร 
NB2 = «Е 5 Le у! КЕ mur + mwp) 
NBS = (Y _Wv + Yiy]y!vlv)/1 
NF = NFl + NF2 + NF3 

where NFl = L ฝ ฝ + (Iy-Ix) pq 
NF2 = - (М мр + ТӨРІМЕ 1 
NF3 = -(N му + NIvlvlvlv /1° 


Again setting terms that include w, p, and q (heave, roll 
and pitch) equal to zero, and set number of known terms 
Peete table II into section 1 of CSMP program (unknown co- 


efficients set equal to zero). 
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TABLE II 


Hvarocynamic Coefficients of Ship "р" for Nonlinear 
Terms (non-dimensional) 


Хо = 0.00005 
Ху = 0.00241 
X - -0.00341 
VV 
Y. |= -0.0416 
“|| 


Rv | v|* -0.01002 


PROPULSION RATIO A η 


ก 0 5 ο ο 1 
А1 -0.00004 -0.00032 -0.00117 
А2 -0.00035 0.00070 -0.00100 
АЗ 0.00099 0 -0.00085 


Fig. 8 - Fig. 12 are the same results as Fig. 3 - Fig. 7. 
The results of computer program 2 are more accurate than 
the computer program 1, when compared with free running 
model test Of NSRDC {Ref. 6]. Fig. 13 and Fig. 14 when 
studying the stability of the ship by applying a force 
Bent to the ship in computer program, set 

NFl 0.0001* (STEP (10.0)-STEP(10.01)) (can use ΝΕῚ be- 
lause NFl in this program equal zero) 

Study direction stability of the ship by plotting 

direction of the ship (used advance VS. transfer ship 
length)and check heading angle of the ship by plotting 


YAW VS. TIME, 
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IV. CONCLUSIONS 


The equation of motion of surface ship and computer 
program developed here including all of six degrees of 
freedom, but the study in III concerns only three degrees 
of freedom (surge, Sway, yaw) because hydrodynamic coeffi- 
cients are not available; when the state of the art 
reaches the stage in which hydrodynamics coefficients are 
available, this computer program can be used in all six 
degrees of freedom. 

Some results from III are not too perfect because the 
lack of some constants and coefficients such as the value 
of mass (m), initial velocity (ADOTO), command speed (UC), 
etc. But for study can adjust from curve for the model 
test [Ref. 6]. 

This computer program did not include some external 
effects such as effects of wave and wind, but these effects 
could be included in the program by adding terms to the IF 
equation. 

The following implementations ME estes for the 
future work. 

A. Study all six degrees of freedom. 

Б; rer the effects of waves and wind. 

C. Study for control of the velocity and direction 

of the ship (by use of "MACROS', "PROCEDURE" or 


subprogram in CSMP). 
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